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Abstract We calculate the coefficients of operators with dimensions d < 7 in the 
operator product expansion of correlators of qTQ currents, for the effective field 
theory of an infinite-mass quark, Q. Exact two-loop results are obtained, with 
an arbitrary gauge group and spacetime dimension, for the perturbative {d = 0) 
and quark-condensate {d = 3) contributions, confirming our previous result for the 
anomalous dimension of the current. Leading-order results are given for light-quark 
operators with d = 5,6,7 and gluon operators with d = 4, 6. The existence of a 
perturbative correction of order 100% precludes a reliable determination of /b from 
non-relativistic sum rules. 
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1 Introduction 



The effective field theory (EFT) of a static quark ||l], ^ considerably simplifies the 
QCD analysis of hadrons containing a single heavy quark. In particular, mesons 
in EFT are analogous to the hydrogen atom in QED, rather than to positronium. 
To leading order in 1/M, EFT possesses a heavy-quark spin symmetry Q, enabling 
hadrons to be classified according to the angular momentum and parity, , of 
light fields. The Qq ground state with j'^ = ^ (S'-wave antiquark) gives 0~ and 
1~ mesons with identical properties (me = m-Q*, /b = /b*)- Excited states with 
jn _ 1+ 3+ (^p_-^ave antiquark) give degenerate pairs of mesons O"*", 1"*" and 
1+, 2+. The splitting in each pair (i.e. the hyperfine structure) is of order 1/M. 

In QCD, however, even such a simplified problem requires nonperturbative 
methods, such as lattice simulation Q. An alternative method is provided by sum 
rules based on the operator product expansion (OPE). 

In order to investigate meson properties in EFT, we consider correlators of 
bilinear currents, J = qTQ, where the tildes distinguish EFT quantities from those of 
conventional QCD and the static-quark field satisfies Q = 'joQ- The irrelevance of 70 
reduces the time-like component of the (pseudo)vector current to the (pseudo) scalar 
current; aoi to 7^; aij to eijfc7fc75- There are thus four such currents in EFT, with 
r = 75, 7j and F = 1, 7^75. The quantum numbers of the first pair, with F 
anticommuting with 70, are those of the ground-state j'" = ^ mesons; the quantum 
numbers of second pair, with F commuting with 70, are those of the excited-state 
jTT _ 1+ j^esons. Currents with quantum numbers of other mesons necessarily 
involve derivatives. 

In the MS scheme, the renormalized current J{fJ-) is obtained from the bare 
current Jq = (JqTQq, in D = 4 — 2e dimensions, by the multiplicative renormalization 
J(/i) = Jl'^^Zj {iJ.)Jo, with jj"^ = ^'^e^ 14:71. The renormalization constant Zj was 
calculated at the one-loop level in Q and to two loops in I?], |8| . The EFT currents 
J are related to the corresponding QCD currents J at the one-loop level by [0 



J(//) = J{p) 

0(«s(/")) 



^ aJu) f H'^ -10 M2 3 , ,1 ^ 

1 + Cf^^ log -^--H^ + HH' T + 4 

47r \ 4 /i2 4 2 J 
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where 7/jF7^ = HT, H' = dH/dD, and in this and all subsequent equations we 
adopt the convention that the upper sign corresponds to F anticommuting with 70 
and the lower sign to the commuting case. Each EFT current, qTQ, is related to two 
QCD currents, qTQ and IjTjqQ, by formula (|l|), which we have verified by methods 
simpler than those in using dimensional regularization of infrared singularities, 
as in [^. 

In co-ordinate space, the correlator i<TJ{x)J^{0)> is proportional to (5(x). Its 
Fourier transform 

n(go) = ijdx e'''^<rj(x)jt(0)> (2) 
therefore depends only on the energy go and is given by the dispersion relation: 



p{uj')duj' 
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where the dots denote a quadratic subtraction polynomial. The contribution of the 
lowest meson, M, is pm(w) = /^'^(w — ujm), where |<0| J(/x)|M>| = /m(/^) is its 
coupling, and ujm its binding energy. Note that the EFT meson state |M> lacks the 
relativistic normalization factor \/2E of QCD and that a; and wm, like all energies 
in EFT, are measured relative to the pole mass M of the heavy quark. 



The correlator (0) has the structure 



n(w) = Tr 



2{^HTi?M}, (3) 



where T = 7or^7o and A and B derive from operators with even and odd dimen- 
sions, respectively. Thus the correlators of ^75(5 and q'jiQ coincide, yielding the same 
ground-state sum rule. Similarly, the correlators of qQ and q^ij^Q coincide, yield- 
ing an excited-state sum rule obtained by a change of sign of the odd-dimensional 
contributions. 

Mesons with a single heavy quark have been considered in applications of nonrel- 



ativistic |ri|, relativistic Borel-transform [l^ , |T4| , and moment-method |15, 



IC, QCD sum rules. (A comparison is made in ||T8|.) EFT sum rules are similar 
to nonrelativistic [^] sum rules, but in the QCD case the perturbative corrections 
are plagued by powers of a large hybrid logarithm: log(Af/a;). In EFT these are 
summed, ab initio, by the introduction of a new anomalous dimension, leading to 
well defined radiative corrections, whose size has so far not received due attention. 
Here we remedy that state of affairs. 

In Sections Q and |^ we calculate the perturbative {d = 0) and quark-condensate 
{d = 3) contributions, up to two loops, using the method proposed in 0. These con- 
tributions should dominate the sum rule. In the Section Q we give higher-dimensional 
contributions, up to d = 7, to leading order. A Borel-transform EFT sum rule 
is derived in Section ^ and conclusions are drawn from a numerical investigation. 
Throughout, we omit 1/M corrections, though these may be systematically included. 
We believe that a proper account of the radiative corrections in the M — > 00 limit is 
required first, in order to determine whether the "QCD hydrogen atom" is amenable 
to sum-rule analysis. Our conclusion has importance for any more realistic analysis 
of the B meson. 



2 Two-loop perturbative contribution 

In we gave a method for calculating two-loop EFT diagrams, with one exter- 



nal momentum and zero light-quark mass, and implemented it as a REDUCE |1S] 



package that evaluates these diagrams in terms of the basic structures ri^2) where 

-2L^)-2^r(-e)/(4^)^/2] " r(l + 2ne). 
Using this package, we evaluate the bare perturbative correlator, up to two loops: 



Tn 



(4) 



For the sum rule, we need the spectral density, which is easily found as the 
discontinuity of (0) divided by 27ri. For the sake of reliability, we have also derived 
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it by direct use of the Cutkosky rules in D dimensions, which yields 



S2 



(5) 



r(l-e) 
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(47r)^/2 
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2C7f5^(2u;)-2- (1 - e)(l - e/2)r(l - 2e)r(l - e) 
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where Pq^^ is the leading-order (one-loop) result, 62 is the (complex) one-loop vertex 
correction to the dia gram with the two-particle cut, and ^3 is the contribution of 
three-particle cuts. Expressing gQ in terms of as(/^) and using Zj we obtain, by 
each method, the renormalized spectral function 



2^ 



1 + 3Cf 



Qs(/^) 
4tt 



, 2uj 47r2 17\ ^ / o, , 
-21og — + — + y j +0[al(ji) 
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which we now compare with the results of more laborious QCD calculations. 



The perturbative QCD (pseudo)scalar correlator, for any pair of quark masses, 
was calculated up to two loops by one of us in |2^]. The (pseudo)vector correlator 
was found in [^], where the results of calculations of absorptive parts |^ were 
compared and corrected. (See also [|l^].) The pseudoscalar and pseudovector results 
were obtained in the scheme with 75 anticommuting with all 7^, by changing the 
sign of a quark mass in scalar and vector results. We may therefore obtain EFT 



results from p3| , 2C] by: transforming the heavy-quark MS mass to the pole mass 
M; using (|^) in the anticommuting 75 scheme; taking the limit M — > 00, with u) and 
m fixed, where = (M -|- u>)'^ and m is the light-quark MS mass. For the spectral 
function this gives 



NcX{uj ± m) 



CFCts{iJ,) J 6a;2 it 6a;m 



18m'' , a 
log — 

m 
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, 2uj 47r2 17 1 ^ / 2/ s 

-41og^ + ^ + 8}+o(a.V) 
-61og^ + 2| + o(a2(^)) +o(m^) 



where Li(2;) 
with (|^), as 



m 



2)^/2_ 'pj^g agreement of 



-Jo dylog{l - y)/y and A 

— > 0, confirms the matching conditions (|l|) for (pseudo) scalar and 
(pseudo)vector currents. The expansion of the full correlator beyond order m? in- 



volves logm terms, which must be absorbed into quark condensates |13|. 
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3 Two-loop <qq> contribution 



One- and two-loop diagrams for the <qq> contribution in any dimension D and in an 
arbitrary covariant gauge can be systematically evaluated by the method developed 
in [^. We have calculated the 11 contributing one-particle- irreducible diagrams 
using our REDUCE package. Combining these with the propagator and vertex 
results of 0, we find a gauge-invariant <qq> contribution to the bare correlator, of 
the form 



ulioo) = ±^^[i + PiriCF5o + {(^AiCA + PFiCF)r? 

+ (i^A2CA + PF2CF + PL2rFiVL)r2}CF5o+0(5o)], (8) 

with the following coefficients, for any gauge group: 

(Z)-l)(Z)-4) 
D-3 

{D^ - 8L>2 + _ 10) (D - 5) 
2(L>-3)2(L>-4) 
2{2D -9){D -2) 
- {D-3){D-A) 

2D'^ - 51D^ + 552D^ - 3245D^ + 11064L>3 - 216261)2 + 22080L> - 8808 
4{2D - 7)iD - 3)2(D - 4){D - 6) 
L>6 _ + mO'^ - 7281)3 + 1704L»2 _ 1968L> + 816 



Pi 
Pai 
Pfi 
Pa2 

Pf2 
Pl2 



2(L'-3)2(D-4)(L>-6) 
4{D - 2){D - A) 



{2D -7){D -3){D -6) 



The infinities in (^) relate the anomalous dimensions of J and <qq>, confirming 
our result ||^ for jj. Using Zj, Zgg, and Z^, we obtain the finite renormalized result 



2uj \ 2 vr 

2 



+ Cf 



vr 




+ OlaM }. (9) 



In QCD, with Ni^ = 4, the ratio of next-next-to-leading to next-to-leading terms 
is large: (— 3.001og(2u;//i) -|- 7.14)Qs(^)/7r. This vanishes for /x = 0.185u;, which is 
too low a renormalization scale to use in EFT sum rules. However, the two- loop 
correction is not large by itself: about 15% for ~ 2a; ~ 1 GeV. We note that two- 
loop corrections to coefficient functions of quark condensates in light-quark QCD 
sum rules were calculated in ||2^ , where the corresponding ratios were found to be 
even larger. 
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4 Higher-dimensional condensates 



We have calculated the contributions of all quark and gluon condensates with d < 
to leading order. 
0, 
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They are most easily found |23] in the fixed-point gauge, 
in which the quark and gluon fields have gauge-covariant Taylor 
expansions and external gluons do not interact with the static quark. Using a sys- 
tematic method |24] to reduce vacuum expectation values to a minimal basis, we 
find the quark-condensate contributions 



<qq> 

2uj 



± 



6<gG 



16u;3 



3<qG 



96u;4 



Gf,ui75q> 



^<qG^\GxvCF^i,q> 



+ 2<qDri,Jui(Tauq> 



(10) 



where G^j. = gG^^J"-, G^^ is its dual, = gJ'^t", and 



DuG'^i^ 



9 Eg qiiit°'q- 

The first three terms agree with |1C] and may also be confirmed by taking the 
M ^ oo limit of the relativistic results in fl^ , |T3| . 

In the fixed-point gauge, no G^ term from the light quark-propagator 5(0, x) 
survives vacuum averaging. Hence there is no one-loop d = 4 gluon-condensate 
contribution in the non-relativistic limit. We obtain the G^ contribution from the 
M — > oo limit of the results of [^] and neglect the term, which is commensurate 
with the unknown radiative correction to the d = 6 quark-condensate contribution 
in (|lO|). The gluon-condensate contribution is thus: 



^ 'i fabcr^ia r^h r~ic ^ 
<5 / '^Xfj.'^iJ.v'^vX^ 



46087r2a;4 



(11) 



5 Sum rule and conclusions 



To obtain a sum rule, we specialize to the SU(3) gauge group of QCD, with Ca = 
A'c = 3, Cp = 4/3, Tp = 1/2, and A^l = 3 or 4 light fiavours. From the two-loop 
anomalous dimension |0] of J and the two-loop /3-function, we have 

/m(/x) = IMa-^/\^^) [l - + O (a,2(M))) , (12) 

5 285 - 77r2 107 , 2^^ 

K= ; ^—TTT, 6=ll--iVL, 

12 276 262 ' 3 

where /m is a renormalization- group invariant. The radiative correction here is 
small: K = 0.189 or 0.227 for Ni^ = 3 or 4. The conventional QCD decay constants 
/m are related by (|I|) to /m, as follows: 

/m = /5/m(M) (l - + O (a.2(M))) + O (^) , (13) 

where C = 4/3 for vector mesons and C = 2/3 for pseudoscalar mesons (if a fully 
anticommuting 75 is used). In analogy with ([T2D, we also have 

<qq>, = <rq>a;'/'{f^){l-K'^ + 0{aUf^))y 
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where <qq> is a (negative) renormalization-group invariant. 



We adopt the standard model of the continuum spectral density, setting it equal 
to the perturbative one, starting at an effective threshold lOq- We use the vacuum 
factorization approximation for d = 6 and d = 7 quark condensates. For the 
latter, this is unreliable but serves as a rough guide, to ensure that the contribution 
is kept small. The d = 5 quark condensate is denoted by mQ<qq>. We omit the 
two- loop term in (^); to use it consistently we would need the unknown three-loop 
term in 7j. Applying a Borel transform ||5|, |l^ to the correlator, we obtain 



f2 p-ujyi/E 



3E^ 



1 



<fg> 
^ 2 



1 + 
1 + 



as(j») 
vr 



UJr 



47r^ 



17 



264 - Utt^ 



m, 



27b 



o_ 

16^2 



+ 



288^4 



TTas<qq>^ 



3 ^ „3 fabcr^a rib ric 



L{x) 



324^3 

\ogx + -i + Ei{x) + i(3 + x)e-^ - I 



(14) 



1 



{l + x + \x'^^ e- 



3 -r ]^gx -r iqqqX -r <^ \x j , 
log X + 7 — I -I- O (e~^x2 log ; 



X <C 1, 
X > 1, 



where the upper (lower) sign is for ground (excited) state mesons, 7 is Euler's 
constant, and Ei{x) = dy e~^/y. 

The left-hand side of sum rule (ID is a renormalization-group invariant; so is the 
right-hand side, to the given order in as{f^). Numerical investigation, over a suitable 
range of the Borel variable E, reveals that all terms may be kept under control, 
save one. The glaring exception is the perturbative correction, which invariably 
exceeds (— 2 log(2ci;c//^) + ll.l)as(jw)/7r. For /i ~ ~ 1 GeV, this is of order 
100%, an order of magnitude greater than in the QCD sum rule for p |5| and, in 
our opinion, far too large for one to have any confidence in the neglect of unknown 
three-loop terms. An alternative way to express this problem is to say that fastest 
apparent convergence, with a vanishing next-to-leading correction, occurs at a scale 
/i < O.OOSwc ~ 10 MeV, an order of magnitude smaller than Aj^. The radiative 
correction to the quark-condensate contribution, on the other hand, is acceptably 
small: around 15%. 

We conclude that non-relativistic sum rules for mesons with a heavy quark are 
intrinsically unreliable. Whilst EFT solves the hybrid log problem, it also leaves one 
with a unique perturbative correction whose numerical size precludes any reliable 
determination of fu and hence any estimate of /b- 
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